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ABSTRACT 

I present simple but robust estimates of the types of sources making up the 
faint, sub-/^Jy radio sky. These include, not surprisingly, star- forming galaxies 
and radio-quiet active galactic nuclei but also two "new" populations, that is 
low radio power ellipticals and dwarf galaxies, the latter likely constituting the 
most numerous component of the radio sky. I then estimate for the first time the 
X-ray, optical, and mid-infrared fluxes these objects are likely to have, which are 
very important for source identification and the synergy between the upcoming 
Square Kilometre Array (SKA) and its various pathfinders with future missions 
in other bands. On large areas of the sky the SKA, and any other radio telescope 
producing surveys down to at least the /zJy level, will go deeper than all currently 
planned (and past) sky surveys, with the possible exception of the optical ones 
from the Panoramic Survey Telescope and Rapid Response System and the Large 
Synoptic Survey Telescope. The Space Infrared telescope for Cosmology and 
Astrophysics, the James Webb Space Telescope (JWST), and in particular the 
Extremely Large Telescopes (ELTs) will be a match to the next generation radio 
telescopes but only on small areas and above ~ 0.1 — 1 //Jy (at 1.4 GHz), while 
even the International X-ray Observatory will only be able to detect a small 
(tiny) fraction of the /xJy (nanoJy) population. On the other hand, most sources 
from currently planned all-sky surveys, with the likely exception of the optical 
ones, will have a radio counterpart within the reach of the SKA. JWST and 
the ELTs might turn out to be the main, or perhaps even the only, facilities 
capable of securing optical counterparts and especially redshifts of /zJy radio 
sources. Because of their sensitivity, the SKA and its pathfinders will have a 
huge impact on a number of topics in extragalactic astronomy including star- 
formation in galaxies and its co-evolution with supermassive black holes, radio- 
loudness and radio-quietness in active galactic nuclei, dwarf galaxies, and the 
main contributors to the radio background. 

Key words: galaxies: active - galaxies: dwarf - galaxies: star formation - radio 
continuum: general - infrared: galaxies - X-rays: galaxies 



1 INTRODUCTION 

The radio bright (> 1 mjy) radio sky consists for the 
most part of active galactic nuclei (AGN) whose radio 
emission is generated from the gravitational potential as- 
sociated with a supermassive black-hole and includes the 
classical extended jet and double lobe radio sources as 
well as compact radio components more directly associ- 
ated with the energy generation and collimation near the 
central engine. Below 1 mjy there is an increasing contri- 
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bution to the radio source population from synchrotron 
emission resulting from relativistic plasma ejected from 
supernovae associated with massive star formation in 
galaxies. After years of intense debate, however, this con- 
tribution appears not to be overwhelming, at least down 
to ~ 50 /^Jy. Deep (Si.4GHz ^ 42 /xJy) radio observa- 
tions of the VLA-Chandra Deep Field South (CDFS), 
complemented by a variety of data at other frequen- 
cies, imply a roughly 50 /50 split between sta r- forming 
galaxies (SFG) and AGN (|Padovani et alj2009lh in broad 
agreement with other rece nt papers (e.g.. ISevmour et al.l 
2008|; ISmolcic et afll200St) . About half of the AGN are 
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radio-quiet, that is of the type normally found in op- 
tically selected samples and characterised by relatively 
low radio-to-optical fl ux density ratios and radio powers 
l|Padovani et alJl2009l ). These objects represent an almost 
negligible minority above 1 mjy. 

This source population issue is strongly related to 
the very broad and complex relationship between star 
formation and AGN in the Universe. At the cosmological 
level, the growth of supermassive black holes in AGN ap- 
pears to be co rrelated with the growth of stellar ma ss in 
galaxies (e.g.. iMerloni. Rudnick fc Pi Matteoll200Sft . At 
the local level, the accreting gas feeding the black hole 
at the centre of the AGN might trigger a starburst. The 
black hole in turn feeds energy back to its surroundings 
through winds and jets, which can compress the gas and 
therefore accelerate star formation but can also blow it 
all away, thereby stopping accretion and star formation 
altogether. Although the details are still not entirely set- 
tled, there is however increasing evidence that in the co- 
evolution of supermassive black holes and galaxies nu- 
clear activity plays a major role through t he so-called 
"AGN Feedback" (e.g.. ICattaneo et alj|2009h . 

Radio observations afford a view of the Universe un- 
affected by the absorption, which plagues most other 
wavelengths, and therefore provide a vital contribution 
to our understanding of this co-evolution. However, while 
we have a reasonably good handle on the radio evolution 
and luminosity functions (LFs) of powerful sources (e.g., 
radio quasars), the situation for the intrinsically fainter, 
and therefore more numerous, radio sources is still murky. 
For example, it is still not entirely clear how strongly 
low-lu minosity, Fanaroff-Riley (FR ) type I radio galaxies 
evolve (|Gendre. Best fc W all 2010). Moreover, the widely 
used SFG LF of I Condon! (119891 ) is highly uncertain at 
the low end (Pi.4GHz < 10 20 W Hz -1 ) because, although 
these sources are intrinsically very numerous, the vol- 
ume within which they can be detected is small (their 
"visibility function" is very low). Indeed, a source with 
-Pi.4GHz ~ 10 20 W Hz -1 could b e detected in the deepes t 
radio image currently available (Owen & Morrison 2008) 
out to z ~ 0.05, but it is not (|Strazzullo et al.ll2010T ). be- 
cause the survey area is too small. And already at z ~ 1, 
assuming a luminosity evolution oc (l + z) 3 , such a source 
would have a flux density as low as ~ 0.2 ^tjy. Finally, 
there are still no published results on the radio evolution 
of radio-quiet AGN, which are intrinsically weak sources 
(fi.4GHz ;$ 10 24 W Hz -1 ) and a non-negligible compo- 
nent of the sub-mjy sky. Note that f aint (sub- /xjy ) radio 
sources have also been proposed by ISingal et all |2010h 
as the main contributors to the extragalactic radio back- 
grou nd recently reporte d by the ARCADE 2 collabora- 
tion |Fixsen et alj|20icl ). Deeper radio observations over 
large areas of the sky are desperately needed to determine 
the LF and evolution of the most common radio sources 
in the Universe. 

These will soon be realised, as radio astronomy is 
at the verge of a revolution, which will usher in an era 
of large area surveys reaching flux density limits well be- 
low current ones. The Square Kilometre Array (SKA j3, in 
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fact, will offer an observing window between 70 MHz and 
10 GHz extending well into the nanojy regime with un- 
precedented versatility. The field of view will be large, up 
to ~ 200 deg 2 below 0.3 GHz and possibly reaching ~ 25 
deg 2 at 1.4 GHz. First science with ~ 10% SKA should 
be near the end of this decade. Location will be in the 
southern hemisphere, either Australia or South Africa. 
Many surveys are being planned with the SKA, possibly 
comprising an "all-sky" 1 /^Jy survey at 1.4 GHz and an 
HI survey out to redshift ~ 1.5, which should include 
~ 10 9 galaxies. 

The SKA will not be the only participant to this 
revolution. The LOw Frequency ARray (LOFAR^Q has 
recently started operations and will carry out large 
area surveys at 15 , 30, 60, 120 and 200 MHz (see 
iMorganti et aU2009l . for details), opening up a whole new 
region of parameter space at low radio frequencies. Many 
other radio telescopes are currently under construction 
in the lead-up to the SKA including the Expanded Very 
Large Array (EVLAjf), the Australian Square Kilometre 
Array Pathfinder (ASKAPfl, the Allen Telescope ArrajQ 
Apertijf], and Meerkafl These projects will survey the 
sky vastly faster than is possible with existing radio tele- 
scopes producing surveys covering large areas of the sky 
down to fainter flux densities than presently available. 

Current deep surveys include a number of VL A small 
area surveys below 0.1 mjy at a few GHz , reaching a 
maxi mum area of ~ 2 deg 2 CVLA-COSMOS; lBondi et alj 
120081 ) and a minimum flux density ~ 15 /iJy at 1.4 GHz 
(SWIRE; [Owen fc Morrisonl|200S ) and ~ 7.5 ^iJy at 8.4 
GHz (SA 13: lFomalont et alfeooa). T he NRAO VLA Sky 
Survey CNVSS; ICondon et al.lll99Sl ) and the Faint Im- 
ages of the Radio Sky at Twenty centimeters (FIRST; 
iBecker et ai1ll995T ) only reach 1.4 GHz limits of 2.5 and 
1 mjy respectively but cover much larger areas (82% and 
22% of the sky respectively). 

What lies beneath the surface of the deepest sur- 
veys we currently have? Predictions for the source pop- 
ulation at radio flux densities < 1 /i J y hav e been made, 
amongst others, by iHopkins et al.l (I2000l) , IWindhorstl 
d2003h, I Jackson! (120041 ), Ijarvis fc Rawlingsl (|2004l ). and 
IWilman et alj (|2008T ). These papers have presented de- 
tailed estimates for the number counts of faint radio 
sources predicting, for example, that SFG should make 
up ~ 90% of the to tal population at Si.4GHz ~ 1 AtJy 
l|Wilman et al. 2008) but had to rely, for obvious rea- 
sons, on extrapolations. Only the last two papers include 
radio-quiet AGN in their modelling by converting their 
X-ray LF to the radio band assuming a linear correlation 
between radio and X-ray powers. Most importantly, as 
described below, two crucial constituents of the sub-/ijy 
sky have been excluded by all of these studies. 

The first aim of this paper is then to have a broad 
look at the likely astrophysical populations, which make 
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up the very faint radio sky. But detecting sources is only 
part of the story, as then comes the identification process. 
This requires a wealth of multi-wavelength data, rang- 
ing from the optical/near-IR imaging needed to provide 
an optical counterpart and, when needed, photometric 
redshifts, to the optical/near-IR spectra required to es- 
timate a redshift, and hence the distance of sources, to 
the X-ray data, which are vita l to separate AGN from 
SFG (e.g.. [Padovani et al.ll2009l . and references therein), 
to the mid-infrared colours, which provide additional in- 
form ation on this separation (e.g.. ISaiina. Lacv fc Scotn 
2005). The second goal of this work is then to provide es- 
timates for the the X-ray, optical, and mid-infrared fluxes 
these sources are likely to have. Optical magnitudes will 
also determine how feasible it will be to obtain redshifts 
for them. This kind of information is important for plan- 
ning purposes, to be ready to take full advantage of the 
new, deep radio data, and also to maximise the synergy 
between the SKA and its pathfinders, and present but 
also, most importantly, future missions. To the best of my 
knowledge, these estimates have never been made before. 

Throughout this paper spectral indices are written 
S v oc v~ a , magnitudes are in the AB system, and the 
values Ho = 70 km s" 1 Mpc" 1 , £l M = 0.3, and tt A = 0.7 
have been used, unless otherwise noted. All the radio re- 
sults refer to 1.4 GHz. Although I mention only the SKA 
throughout the paper, my results obviously apply to any 
other radio telescope producing surveys reaching the fijy 
level. I also make reference to a large number of missions 
and projects, for which I give mainly sensitivity infor- 
mation. Readers wanting to know more should consult 
the relevant World Wide Web pages, whose addresses are 
provided in the text. 



2 MICROJY AND NANOJY RADIO 
SOURCE POPULATION 

I present here a simple approach to study the radio sky 
source population, based on only two parameters: the 
smallest flux density and the largest surface density of 
radio sources. The main idea is to provide robust results 
based on some basic observables and to pay particular 
attention to all populations reaching below the fijy level. 

The smallest flux density /n m of a population of 
sources depends on the minimum radio power at z ~ 0, 
■Pmin(O), the maximum redshift of the sources, z max , and 
any luminosity evolution le(z), where P(z) = P(0) x 
le(z). If evolution peaks at z top z max ) and then stops, 
then 



/lim = 



Pmin (0)Ze(Ztop)(l + Zmax) 1 " 

47rP>|(z max ) 



(1) 



where Dl(z) is the luminosity distance. Luminosity evo- 
lution makes sources brighter and therefore increases /ii m - 
The largest surface density of a population, N(> 
/lim), depends on its number density at z ~ 0, iVr(0), 
the maximum redshift z max , and any density evolution 
de(z), where N T (z) = Nt (0) x de(z). Then 



where dV/dz is the derivative of the comoving vol- 
ume. In case of no density evolution N(> /n m ) = 
ATt(0)V(w)/47t sr" 1 . 

Information on the local LF, needed to derive Nt{0) 
and Pmin(0), and the evolution for various classes was 
derived from a variety of sources. Namely: 



(i) the LF for SFG is that of ISadler et"aH (I2002TI 
which at the low end is complemented by that of lCondonl 
(1989); luminosity and density evolution and z top are 
from lHopkinsl i|2004l ); 

(ii) the LF for ra d io-qu iet AGN is built using 
data from iRush et al.1 (Il996l ) (12/jm and CfA Seyfert 
samples) and Padovani et al. (in preparation: VLA- 
CDFS sample). Evolution is from Padovani et al. 

in preparation), while Ztop is fr om the X-ray band 
Hasinger. Mivaii. fc Schmidtll2005l l: 

(iii) the LF for F R Is derives from that of 
lUrrv fc Padovanil (1 19951), which agrees well with t he very 
recent derivation of Gendre. Best fc Walll ((2010) apart 
from the first bin; the LF was then modified accord- 
ingly and the values of Nt( 0) and P m in(0) were derived. 
iGendre, Best fc Walll \20Vj ) find evidence of evolution for 
FR Is but only at Pxaghz > 10 25 W Hz -1 . Since the to- 
tal surface density in eq. [2] depends on the most numer- 
ous, and therefore least luminous, sources (which reach 
Pi.4GHz ~ 10 23 W Hz^ 1 ) no evolution was assumed; 

(iv) the LF fo r FR lis is also based on that of 
lUrrv fc Padovanil (Il995h , whi ch ag rees quite well with 
that of lGendre. Best fc Walll (|2010h down to Pi. 4 ghz ~ 



10 W Hz . Below this value the LF is basically un- 
determined so no modification was done. Evolution was 
taken from lUrrv fc Padovanil (Il995l ); 

(v) the LF for flat-spectrum radio quasars (FSRQs), 
steep-spectrum radio quasars (SSRQs), and B L Lacs 
are based on those of lUrrv fc Padovanil (|l995l ) , which 
were calculated from those of FR lis and FR Is re- 
spectively, based on a beaming model. These have been 
shown to agree with LFs derived from recent samples 
(|Padovani et al.ll2007h . Since the FR I LF was modified, 
the beamed LF for BL Lacs was also changed accordingly 
at low powers before deriving Nt(0) and P m in(0). Evo - 
lutionary parameters are from lUrrv fc Padovanil jl995h . 
while ztop values fo r FSR Qs, SSRQs, and FR Ilfl come 
from lde Zotti et al.1 (120051 ). 



Values of z max were fixed to the highest redshift 
of the class under consideration, that is ~ 6, 6.5, and 
5.5 for SFG, radio-quiet AGN, and FSRQs, SSRQs, and 
FR lis respectively, apart from FR Is and BL Lacs, for 
which z max = 3 was a s sume d. Finally, the blazar cata- 
logue o f Massar o et al.l (120091 ) and the AGN catalogue of 
IPadovani et al l l| 19971 ) were checked to see if any sources 
had radio power below the adopted P m i n values. None 
was found. 

Table [T] summarises the parameters used in eqs. [T] 
and[51 JVt(0) values were derived in most cases from sim- 
ple fits to the LF and should be considered approximate. 



JVt(0) 
4-k 



de(z)dV/dz 



(2) 



8 Ac cording to unified schemes (see, e.g., lUrrv fc Padovani 
1995), SSRQs and FR, lis, being FSRQs seen at larger angles 
with respect to the line of sight, need to share the same ztop. 
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Figure 1. The largest surface density vs. the smallest flux density for various classes of radio sources. The two horizontal lines 
denote, from top to bottom, the surface density of the optical sources in the Hubble Ultra Deep Field and the surface density of 
the AGN needed to explain the X-ray background. See text for more details. 



Table 1. Radio populations parameters 



Class 


JV T (0) 
Gpc" 3 


■Pmin(O) 

W Hz" 1 


LE 


DE 


%top 


-2 max 


FSRQs 


12 


2 x 10 24 


exp[T(z)/0.23] a 




2.25 


5.5 


SSRQs 


59 


3 x 10 24 


exp[T(z)/0.15] a 




2.25 


5.5 


FR lis 


590 


3 x 10 24 


exp[T(z)/0.26] a 




2.25 


5.5 


BL Lacs 


2,310 


10 23 


exp[T(z)/0.32] a 






3.0 


FR Is 


29,300 


10 23 








3.0 


RQ AGN 


3.9 x 10 5 


5 x 10 19 


(1 + *) 2 ' 4 




1.7 


6.5 


SFGs 


4.5 x 10 7 


2 x 10 18 


(l + z) 2 - 7 


(! + * a )0.16 


2.0 


6.0 


Dwarf Galaxies 


2.0 x 10 8 


< 2 x 10 18 


(1 + z) 2 - 7 




2.0 


3.0 


Low-power Ellipticals 


4.8 x 10 6 


< 3 x 10 19 




(1+^)" 1 - 7 




3.0 



Hq = 50, qo = 0; T(z) is the look-back time 
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The N t(0) and P m i n (0) values from lUrrv fe PadovanH 
(1995) have been converted to Ho = 70 from Ho = 50 
while the evolutionary parameters still refer to an Ho = 
50, go = cosmology (as a simple conversion in this case 
is not possible). Note that, although more complex evolu- 
tionary models than those used here have been adopted 
in the literature for "classical" powerful radio sources, 
the consen sus is that such sources reach only the ~ mjy 
level (e.g., |jacksonll2004l ; IWilman et al.ll200Sl V This is in 
agreement with my own results below and in any case 
well above the flux densities of interest here. 

The resulting flux and surface density limits are 
shown in Fig. [1] These values should be considered as 
robust upper and lower limits respectively, because: 1. 
one cannot exclude that lower-power, and therefore more 
numerous, objects exist; 2. z max could be larger than as- 
sumed. The latter parameter has a stronger influence on 
flux than on surface density. For example, for z max = 6, 
eq. [T] shows that the former value decreases by a fac- 
tor (4/7) 1 ' ar [D L (6)/D L (3)] 2 ~ 4.4 (a r = 0.7) as com- 
pared to z max = 3, while from eq. [2] and in the case 
of no density evolution the latter increases by a factor 
V(6)/V(3) ~ 2.2. If z m ax increases from 6 to 10, the two 
values are instead ~ 2.8 and ~ 1.5. 

Fig. Q] shows that the most powerful radio sources, 
that is FRSQs, SSRQs, and FR lis are, not surprisingly, 
the ones having the largest flux density (ss 0.1 — 1 mjy) 
and the smallest surface density (~ 1 — 50 deg" 2 ) lim- 
its. BL Lacs are only slightly fainter then FR lis, while 
FR Is are the only radio-loud sources reaching « 1 fijy. 
Radio-quiet AGN and SFG are the faintest classes, going 
into the nanojy regime, with SFG dominating the faint 
radio sky (amongst "classical" rad io sources: see b elow) . 
Indeed, the differential counts of IWilman" et all (|2008h 
predict a strong dominance (~ 90%) of SFG at 1 fijy, 
which gets slightly weaker at nanojy flux densities. 

In Fig. [T] I have also plotted two horizontal lines. 
The one at N w 10 4 deg -2 denotes the surface density of 
the AGN needed to explain the X-ray background (Gilli, 
private communication, based on lGilli et al.l l2007). Since 
these for the most part are radio-quiet AGN, its proxim- 
ity to the estimated surface density of this class is reas- 
suring. The other line at N ~ 2 x 10 6 deg -2 indicates the 
surface density of t he optical sources in th e Hubble Ultra 
Deep Field (UDF) l|Beckwith et alj|2006h . thought to be 
mostly star- forming galaxies. The fact that this is very 
close to the completely independent estimate for SFG in 
the radio band shows that the latter is perfectly plausible. 

I argue that two other populations play a ma- 
jor role at 5i.4GHz < 1 ^Jy- The first one is that 
of low-power ellipticals. It has been know for quite 
some time that ellipticals of similar optical luminos- 
ity vary widely in radio power, with some (non-dwarf) 
galaxies having Psghz < 2 x 10 19 W Hz -1 (e.g., 
ISadler. Jenkins fc Kotanvl 1 19891. convertin g from their 
cosmology). More recently. ICapetti et alj (|2009h have 
shown that 82% of early-type galaxies in the Virgo clus- 
ter with Bt < 14.4 are undetected at a flux density limit 
of ~ 0.1 Jy, which implies core radio powers Pg.4GHz < 
4 x 10 18 W Hz" 1 . iMiller et all |2009T l have studied the 
radio LF in the Coma cluster and found that 58% of 
red sequence galaxies with M r —20.5 are undetected 



at about 28 fiJy r.m.s. Stacking these sources, they ob- 
tained a detection corresponding to Pi.4GHz ~ 3 x 10 19 W 
Hz -1 . These faint ellipticals are not represented in pre- 
vious models of the sub-/xjy sky: f or example, th e lower 
limit of the radio-loud AGN LF in lWilman et al.l (I2008T ) 
is Pi.4GHz = 2 x 10 20 W Hz" 1 . We have no informa- 
tion on the evolution of these radio sources but, based on 
the results of low-power (Pi.4GHz < 10 25 W Hz" 1 ) radio 
galaxies (e.g.. iGendre. Best fc Waill2010l . and references 
therein), it seems plausible to assume no luminosity evo- 
lution. Taking P min < 3 x 10 19 W Hz" 1 and z max = 3 
one derives fu m < 0.6 nanojy (assuming the same lumi- 
nosity evolution as found by ICimatti. Daddi fc Renzirn] 
(|2006P l for ellipticals in the B-band one would instead get 
/iim < 5 nanojy). (Note that since we only have an up- 
per limit on P m i n in this case /ii m is an even more robust 
upper limit than for the previously discussed classes.) 

As regards their surface density, I have taken the 
n umber density of all early - type galaxies from the LF 
in lde Lapparent et al.l l|2003h (the two- wing Gaussian fit 
for R c ^ 21.5 in their Tab. 7), which agrees within 
~ 20% with the Schechter fit to th e local Slo a n Dig - 
ital Sky Su r vey (S PSS) LF done by iBell et all (|2004 ). 
iKriek et al.l {2008) have recently suggested a decrease 
of a factor ~ 8 in the number density of high-mass 
(> 10 1:l Mq) early-type ga l axies between z ~ and 2.3. 
ICimatti. Daddi fc Renzirn] (|2006T l have found a similar 
trend for low-mass (< 10 n MQ) early-type galaxies be- 
tween z ~ and 1.2. Assuming that all ellipticals are 
radio sources at some (very low) level, such a density 
evolution, and z max = 3, I get from eq. [5] a limiting sur- 
face density for low-power ellipticals « 2.4 x 10 4 deg" 2 , 
of the same order as that of radio-quiet AGN. 

The other population missing from previous studies 
is that of dwarf galaxies, which are very faint and con- 
stitute the most numerous extragalactic population. This 
class includes dwarf spheroidals and ellipticals, dwarf ir- 
regulars, and blue compact dwarf galaxies (BCDs), and 
it has never been considered for the simple reason that its 
radio LF and evolution has never been determined. But 
the simple approach adopted here can provide us with 
some idea of how faint and how numerous these sources 
are going to be in the sub-^tjy sky. 

The determination of the LF of dwarf galaxies is 
not an easy task as it requires a thorough understand- 
ing of selection effects, due to their low surface bright- 
ness. The LF appears also to be environment-dependent 
(e 

LF of 

brightness incompleteness, and derived the number den- 
sity of dwarf galaxies (defined by M r < —18.8 [equivalent 
to Mb < —18.1], which is where there is an upturn in 
the slope of the LF). Available data are consistent with 
no density evolutio n at the faint end of the LF up t o 
at least z w 3 (e.g.. ICooray||2005l ; ISalimbeni et al.ll2008h . 
With these assumptions I derive from eq. [2] a (likely) con- 
servative limiting surface density ~ 5 x 10 6 deg" 2 , higher 

than all other classes. 

As regards flux density, iLerov et alj l|2005l ) have 
shown that P m in,i.4GHz for dwarf galaxies is < 1.6 x 10 18 
W Hz" 1 . Sin ce most galaxies at the faint end of the L F 
are blue (e.g jBfanton et al.|[2005l : I Salimbeni et aJj|2008h . 



Ferguson fc Sandagej|199ll ). I have taken the SDSS 
Blanton et all (12005). which is corrected for surface 
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most dwarfs in the Universe should be of the star-forming 
type and I then assume the same luminosity evolution as 
for SFG. For z max = 31 then get /i im < 0.6 nanoJy (no 
evolution would imply /n m < 0.03 nanoJy). (As was the 
case for low-power ellipticals, since we only have an upper 
limit on P m i n /ii m is an extremely robust upper limit.) In 
summary, dwarf galaxies are likely to be the most numer- 
ous component of the faint radio sky. 

Although this paper deals with extragalactic sources, 
one might worry that stellar objects could also contribute 
substantially to the /j,3y and nanoJy sky. This is ex- 
tremely unlikely. The radio thermal component of the 
Sun at the distance of a Centau ri would have a flux den- 
sity ~ 5 — 30 /iJy ((White 2004) (where the lower value 
refers to the quiet state and the higher one to the flaring 
one), which means ~ 0.7 — 4 nanoJy at 10 Kpc. More- 
over, a Sun-like star at 10 Kp c would have a non-thermal 
flux density ~ 0.001 nanoJy (|Seaquistll 19971 ). Given that 
the most common main-sequence stars are of the M type, 
which are more than one order of magnitude less lumi- 
nous than the Sun, the bulk of stellar radio emitters will 
be very faint (< 1 nanoJy). 



3 MULTI- WAVELENGTH PROPERTIES OF 
MICROJY AND NANOJY RADIO 
SOURCES 

I estimate here the X-ray, optical, and mid-infrared fluxes 
radio-quiet AGN, SFG, and FR Is should have at the /xjy 
and nanoJy flux density levels. Low-power ellipticals and 
dwa rf galaxies require some d iscus sion. 

iBalmaverde fc Capettil (|2006l ) have studied the 
multi-wavelength characteristics of very low radio power 
ellipticals with 10 19 < P 5G Hz < 3 x 10 24 W Hz" 1 
((Psghz) ~6x 10 21 W Hz" 1 ). These sources can be 
considered as miniature radio-galaxies, in the sense that 
their nuclear properties are scaled down versions of those 
of low-luminosity, FR I radio galaxies. Indeed, if one sep- 
arates radio galaxies on the basis of their nuclear ac- 
tivity into high-excitation (HERGs) and low-excitation 
(LERGs) radio-galaxies, almost all FR Is are LERGs and 
most FR lis are HERGs, alth ough there is a po pulation 
of FR II LERGs as well (e.g., lLaing et al.1ll994T l. In this 
scheme, low-power ellipticals could be considered as the 
natural extension of LERGs to lower radio luminosities. 
In the following I will then assume that low-power el- 
lipticals have the same multi-wavelength properties as 
their higher-power relatives, with the obvious caveat that, 
for the same galaxy optical magnitude they will have a 
much lower radio emission, and therefore will be char- 
acterised by a much lower radio-to-optical flux density 
ratio. The situation for dwarf galaxies is more complex. 
In dwarf spheroidals and ellipticals very little, if any, 
star formation is going on now, alt hough their star for- 
matio n histories are complex (e.g., iFerguson fc Bin ggcli 
1 19941 ). Dwarf irregu lars appear to be low mass v ersions 
of large spirals re.g.. lKlein|[l98r3 : lLerov et al.ll2005l ). while 
BCDs fall, on average, at th e high star formation rate end 
ijHunter fc Elmegreer]|2004l ). 

Since dwarf galaxies are intrinsically weak, our 
knowledge of their spectral energy distributions (SEDs) is 



quite scanty and based on small samples, which are most 
likely affected by selection effects. My working hypothe- 
sis, which is not contradicted by the data (see below) , will 
be that dwarf spheroidals and ellipticals have SEDs simi- 
lar to those of low-power ellipticals, while dwarf irregulars 
and BCDs are mini-spirals. As mentioned above, the faint 
radio sky should be dominated by the star-forming type 
of dwarfs. 

In the following the multi-wavelength fluxes of the 
various classes of faint radio sources are estimated us- 
ing three (two for the near-IR band) different methods. 
Readers not interested in the details should skip directly 
to Sect. 13X31 [3331 and [3X21 



3.1 X-ray Band 

3.1.1 Typical X-ray-to-radio flux density ratios 

iRanalli et al.l (|2003T ) have shown that X-ray and radio 
powers in SFG are strongly correlated, likely because 
they both trace the star formation rate. Their typical 
X-ray-to-radio flux ratio is (log(/ . 5 -_2kcv/Si.4Giiz)) = 
— 17.9 ± 0.3 (where the X-ray flux is in c.g.s. units and 
the radio flux density is in ^iJy). Although derived for 
local galaxies, the co rrelation appears to hold at higher 
redsh ifts as well fe.g.. lRanalli et al]|2003l ; IPadovani et al.l 
120091 ). 

lOtt. Walter fc Brinks! (120051 ) have studied the 
X-ray properties of eight dwarf galaxies under- 
going starburst observed with Chandra. I derive 
(log(/o. 5 -2kcv / 'Siaghz)) ~ -18 (where I have used ob- 
served X-ray fluxes converted to the 0.5 — 2 keV band, 
taking into account both extended and point source emis- 
sion, and radio data from the literature), perfectly con- 
sistent with the mean value for SFG. 

X-ray and radio powers are also correlated in radio- 
quiet AGN, with a relatively large dispersion (see, e.g., 
Fig. 13 of iBrinkmann et all |2000| ). In deriving typi- 
cal X-ray-to-radio flux density ratios for this class, it 
is important to kee p in mind that, as discussed in 
I Padovani et al.l |2009l ). X-ray selection will tend to favour 
sources with relatively large X-ray-to-radio flux ratios 
while the opposite will be true for radio selection. I 
have use d the VLA-CDFS rad io-quiet AGN sample se- 
lected bv lPadovani et ah! (|2009l ) (and refined by Padovani 
et al. in prep a ration ) and the X-ray fluxes provided 
by iTozzi et all (|2009l ) to obtain a K-corrected value of 
(log(/o. 5 -2kcv/Si.4GH 2 )) = -17.2±0.9for radio-selected, 
radio-quiet AGN. (For a x ~ 0.9 — 1.1, the range cov- 
ered by (unabsorbed) RQ AGN and spirals, and syn- 
chrotron emission (a r ~ 0.7), K-correction effects al- 
most cancel out.) Note that these radio-quiet AGN in- 
clude both broad-lined (type 1) and narrow-lined (type 
2) AGN. Whenever possible, X-ray fl uxes were corrected 
for absorption fsee lTozzi et al.ll2009l . for details). As re- 
g ards X-ray selected, ra dio quiet AGN, I use the results 
of lPadovani et all J2009h on the hard X-ray selected sam- 
ple of lPolletta et alj (|2007l ) (/ 2 _iokeV > KT 14 erg cm" 2 



s _ ) derived using survival analysis due to the many up- 
per limits on radio flux densities. By converting them to 
the 0.5 — 2 keV band I obtain a K-corrected mean value 

(log(/ .5-2kev/Sl.4GH*)) ~ -14.8. 
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X-ray and radio emission correlate in radio galax- 
ies as well, a l beit w ith a large scatter. Fig. 5 of 
iPadovani etall l|2009l ) shows that radio galaxies have 
an X-ray-to-radio power ratio of the same order of that 
of SFG but smaller than that of quasars. I Evans et al.l 
(120061 ) have studied the X-ray properties of the cores 
of 22 low-redshift (z < 0.1) radio galaxies and find 
a strong correlation between 1 keV X-ray and 5 GHz 
radio core powers for low-absorption (Ah < 5 x 10 22 
cm -2 ) radio galaxies, a sub-sample which includes most 
FR Is. For a x = 0.9 (their average value) and ct r ~ 
(appropriate for radio cores) I get a mean value 
(log(/o.5-2kev/#i.4GH*)) ~ -18. One might worry that 
for more distant sources it will be hard to resolve the 
nuclear components. I then derived the typical X-ray-to- 
radio fl ux density ratio for th e FR Is in the AGN cata- 
logue of IPadovani et all fl997) by using total X-ray and 
radio fluxes, obtaining (log(/o. 5 -2kcv/Si.4Gffz)) ~ -19. 
In the following I will then adopt the intermediate value 
(log(/o. 5 -2kov/S'i.4GH Z )) ~ -18.5, noting that even as- 
suming the largest value would not affect any of my con- 
clusions. 

Flux density ratios will be largely unaffected by evo- 
lution, since this is broadly similar in the X-ray an d ra- 
dio band f or SFG (IRanalli. Comastri fc SettH 120051 ) and 
AGN (e.g.. lWall et al.ll2005h . No information is available 
on the X-ray evolution of FR Is, which in any case are 
weak X-ray emitters. 

This approach has its limitations. The SFG used by 
IRanalli et all d2003l) have 5i. 4G Hz > O.Uy and the FR Is 
studied bv lEvans et all |2006) have S 5 GHz(core) > 0.01 
Jy. The derived mean X-ray-to-radio flux density ratios 
are then used to predict the X-ray fluxes of radio sources 
in the pJy and nanoJy regime, that is at > 10 4 times 
fainter flux density levels. Moreover, these values are 
likely to be prone to selection effects, which I have tried 
to take into account in the case of radio-quiet AGN. The 
latter could be less relevant if these values are based on 
a strong, physically based correlation, which is likely the 
case for SFG. 



3.1.2 Extrapolation from the VLA-CDFS Sample 

An alternative, complementary procedure is to derive the 
mean X-ray flux for a sample with a relatively faint radio 
flux density limit and then shift it down to simulate a 
fainter sample . To this aim, I used t he VLA-CDFS sam- 
ple for which IPadovani et al.l (|2009h (see also Padovani 
et al., in preparation) provide a classification in SFG and 
AGN (radio-quiet and radio-loud). Based on their results, 
most AGN of the radio-loud type are expected to be low- 
luminosity radio-galaxies, that is FR Is. The flux density 
limit of the VLA-CDFS survey is not constant but in- 
creases with distance f rom the field centre from a min- 
imum value of 42 ^ijy l|Kellermann et aLll200St ). It then 
follows that low-flux density sources, which, on average, 
have also the smallest X-ray fluxes, are underrepresented 
as compared to a sample with a constant limit across the 
field. To correct for this I have weighted each source by 
the inverse of the area associated to its flux density. I 
then divided the mean values of the X-ray flux by 42 to 



simulate a sample with Si.4GHz ^ 1 /^Jy. Since most SFG 
and radio-loud AGN are undetected in the X-ray band, 
only upper limits are available for these two classes. In 
an hypothetical sample with a radio flux density limit 
of 1 ^Jy SFG should have (/o.5-2k c v) < 5 x 10~ 18 erg 
cm -2 s _1 , FR Is should have (/o.5-2k c v) < 10~ 17 erg 
cm -2 s _1 and radio-quiet AGN should be characterised 
by (/o.5-2kov) ~ 2 x 10~ 17 erg cm~ 2 s" 1 . 

3.1.3 X-ray fluxes of faint radio sources 

Figure [2] plots 0.5 - 2 keV flux vs. 1.4 GHz radio flux 
density and shows the loci of X-ray selected and radio- 
selected, radio-quiet AGN, SFG, and FR Is (Sect. EXT). 
Note that radio-quiet AGN will span the full range be- 
tween the two dashed lines in Fig. [2] with radio (X-ray) 
selection favouring sources with low (high) X-ray-to-radio 
flux density ratios. The position of these loci with re- 
spect to survey limits determines the fraction of sources 
of a given class detected in one band with counterparts 
in the other. The figure shows also the expected X-ray 
fluxes for "typical" radio-quiet AGN, SFG, and FR Is 
(from the NASA/IPAC Extragalactic Database [NED]) 
scaled to 1 /xjy. The mean radio and X-ray flux values, 
or upper limits, for sources belonging to an hypothetical 
sample characterised by 5i.4GHz ^ 1 £tJy (Sect. ETl~2"j) are 
also shown. The fact that the three estimates give consis- 
tent results is reassuring and shows that we can predict 
reasonably well the X-ray fluxes of faint radio sources. 

The horizontal dot-dashed line at /o.5-2kcV ~ 2 x 
10 -17 erg cm -2 s _1 represents the deepest X-ray data 
currently avai lable, that is the Chandra Deep Field South 
2 Ms Survey i|Luo et al.l 120081 ). covering about 0.1 deg 2 . 
The other horizontal lines indicate, from top to bottom, 
the limiting fluxes for point sources for surveys to be car- 
ried out with eRosit43 and the Wide Field X-ray Tele- 
scopS (WFXT0 The faintest X-ray flux limit corre- 
sponds to the International X-ray Observatory (IXOj3i 
which will provide the deepest X-ray view on the Universe 
for quite some time. IXO will be an observatory type mis- 
sion, like Chandra, and therefore will not produce large 
area surveys. 

The main message of Fig. [2] is that even the most 
powerful X-ray missions we are going to have for the next 
20 years or so will only detect the counterparts of radio- 
quiet AGN with radio flux densities down to ~ 1 /iJy. 
The bulk of the /xJy population, which is most likely to 
be made up of SFG, will have X-ray fluxes beyond even 
the reach of IXO. Same, or possibly even worse, story for 
FR Is. The situation will obviously be even more critical 
in the nanoJy regime, where very few radio sources will 
have an X-ray counterpart in the foreseeable future. 

On the positive side, basically all extragalactic 
sources in the eRosita All-Sky and Wide Surveys and 

9 http:/ /www. mpe.mpg.de/heg/www/Projects/EROSITA/- 
main.html 

10 http://wfxt.pha.jhu.edu/ 

11 A WFXT-ce ntric version of Fig. [2] has been presented by 
iPadovanil 1I2OIOI ) 

12 http://ixo.gsfc.nasa.gov/ 
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Figure 2. 0.5 - 2 keV X-ray flux vs. the 1.4 GHz radio flux density for faint radio sources. The loci of X-ray selected and radio- 
selected, radio-quiet AGN (long-dashed lines), SFG (dotted line), and FR Is (short-dashed line) are indicated. The scaled X-ray 
fluxes of prototypical representatives of the three classes at Si.4GHz = 1 A*Jy are a l so shown. Finally, the mean radio and X-ray 
flux values, or upper limits, for sources belonging to an hypothetical sample characterised by Si.4GHz ^ 1 A»Jy, as extrapolated 
from the VLA-CDFS sample, with error bars indicating the standard deviation, are also marked. The horizontal dot-dashed lines 
indicate the limits of (from top to bottom): the eRosita All-Sky and Wide Surveys, the WFXT Wide and Deep Surveys, Chandra's 
deepest surveys, and IXO. Survey areas and launch dates for future missions, or best guesses at the time of writing, are also shown. 
See text for more details. 



WFXT Wide Survey will have an SKA counterpart, as 
they will have Si^ghz > 1 A*Jy- This should help in the 
identification work of, for example, the 10 million or so 
point sources expected in the latter, by also providing 
very accurate positions. Radio detection of the bulk of 
the AGN in the WFXT Deep Survey will require much 
higher (~ 20 nanojy) sensitivities. This might be accom- 
plished by the SKA given also the small area of the survey 
(~ 100 deg 2 ), under the obvious condition that WFXT 
surveys are carried out in the southern sky. 



3.2 Optical Band 

3.2.1 Typical radio-to-optical flux density ratios 

At variance with the X-ray case, there is in general no 
apparent correlation between radio and optical powers in 
extragalactic sources. This is due to the fact that while 
the radio, far-infrared, and X-ray bands all trace the star 
form ation rate (SFR) (e.g.. lKennicuttl l998: Ra nalli et al.l 
2003), the optical band does not, as young stars domi- 
nate the ultraviolet continuum. Moreover, as mentioned 
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in Sect. [2] elliptical galaxies of similar optical power can 
host radio sources differing by huge amounts in their ra- 
dio luminosity. 

However, SFG and radio-quiet AGN can only 
reach a reasonably well-defined value of the K- 
corrected radio-to-optical flux density ratio R — 
log(5 l i.4GHz/ 5 l R T ^„„), where 5 , r^,„ is the R-band flux 
densi ty (e.g. . [Machalsk r&: Condonlll999l : IPadovani et al.l 
2009). Converting from the values deri v ed fo r the two 
classes in the V-band bv lPadovani et all (|2009h and tak- 
ing the average one gets for the R-band a maximum value 
R « 1.4. Due to K-correction effects (see below), observed 
R values for SFG and radio-quiet AGN can be > 1.4. 
It has to be noticed that, while all "classical" radio- loud 
quasars have R > 1.4, this is not th e case for many radio- 
galaxies (e.g.. IPadovani et alJ l2009). which can extend to 
R < 1.4. Indeed, as mentioned above, low-power ellipti- 
cals will have lower R values than FR Is. 

The absolute SFR in star-forming galaxies spans a 



ies and up to ~ 1000 Mp, yr -1 in the most luminous 
IR starbursts (|Kennicuttl 1 19981 ) . Therefore, for the same 
optical magnitude "normal" spirals will be characterised 
by lower radio emission than starbursts, and therefore 
will have a smaller radio-to-optical flux density ratio. 
To better quantify this I will define as a "normal", non- 
starburst source a galaxy with an SFR < 10 Mq yr -1 , 
which is the maximum value reache d by local Uppsal a 
Galaxy Catalogue (UGC) galaxies (|james et al.ll2004l ). 
Conver ting this to a radio p o wer f ollowing the calibra- 
tion of ISargsvan fc Wedmanl l)2009l ) , which agrees very 
well with that based on polycyclic aromatic hydrocar- 
bons (PAH), translates to Pxaghz < 8.4 x 10 21 W Hz -1 . 
SFG in the VLA-CDF S sample below this value have 
R < 0.3 (and z < 0.2). IPadovani etafl (|2009h have also 
shown that R decreases with decreasing Pi.4GHz for the 
spirals and irr egulars in their sample. From Fig. 1 of 
I Condon! (1 19891 ) (converting from his cosmology) I derive 
R < 0.1 for Pi.4GHz ~ 8 x 10 21 W Hz" 1 for optically 
bright (Bt ^ 12) spirals and irregulars; lower radio pow- 
ers have smaller R values. 

Dwarf galaxies have also low flux density ratios 
and are in fact unde rrepresented in radio samples 
JVan Duvne ct al. 2004). From the mean values given by 
iLerov et al.l ( 2005 ) I get R sa —0.6 ( conve rting to my 
notation), while Fig. 9 of iLeon et alj l|2008h shows that 
R < -0.8 for Mb ~ -17. Even BCDs, which are the 
most star- forming amongst dwarfs, are cha racterised by 
R » -0.2 (|Hunt. Bianchi fc Maiolinoll2005l ). 

In the optical band the K-correction is quite im- 
portant. For example, at z ~ 1, the mean redshift of 
the VLA-CDFS sample, it reaches ~ 1 magnitude for 
Sbc galaxies and ~ 2 magnitudes for ellipticals (e.g., 
IColeman et alJ ll98C0. These values should be compared 
to the factor ~ 1.2 in flux (equivalent to ~ 0.2 mag- 
nitudes) expected in the radio band (for a r ~ 0.7). 
Moreover, absorption by the intergalactic medium will 
also further decrease the optical flux. It then follows 
that the derived magnitudes are very robust lower lim- 
its. As was the case for the X-ray band, flux den- 
sity ratios will be largely unaffected by evolution, since 



this is b roadly similar in the optical an d radio band 
for S FG (iRanalli Comastri fc Sett J [20051 ) and powerful 
AGN (|Wall et al.ll2005l ). For FR Is and low-power ellipti- 
cals, however, the evolutiona ry correction in the optical 
band can be quite important (Poggianti 1997). This, cou- 
pled to the almost absent radio evolution (Sect. [2]), could 
imply a decrease in radio-to-optical flux density ratios for 
these classes at high redshifts. 



3.2.2 Extrapolation from the VLA-CDFS Sample 

The best estimate of the optical magnitudes of faint radio 
sources is obtained by scaling those of a sample with a rel- 
atively faint radio flux density limit. Using the same pro- 
cedure as for the X-ray band (Sect. l3~1.2|) . I added to the 
(-Rmag) of the VLA-CDFS sources 2.5 x log 42 to simulate 
a sample with Si.4GHz ^ 1 A*Jy. Correcting for the non- 
uniform flux density limit as done before (Sect. 13.1.21) . 
SFG in this hypothetical sample should have (R m a g ) ~ 
26.3, radio-quiet AGN should have (-Rmag) ~ 28.2 and 
radio-loud AGN (mostly low-luminosity radio-galaxies) 
should have (-Rmag) ~ 27.3. Note that ~ 85% of the SFG 
in the VLA-CDFS sample have radio powers above the 
"normal" spiral limit and so are mostly of the starburst 
type. 



3.2.3 Optical magnitudes of faint radio sources 

Figure [3] plots R mag vs. 1.4 GHz radio flux density and 
shows: 1. the expected R m a g for "typical" radio-quiet 
AGN, SFG (starbursts, spirals, and dwarfs), and FR Is 
(from NED) scaled to 1 /xJy; 2. the maximum value for 
SFG and the approximate dividing line between radio- 
loud and radio-quiet AGN (diagonal dashed line; see Sect. 
I3.2.1[) . SFG and radio-quiet AGN are expected to popu- 
late the top left part of the diagram; 3. the mean radio 
and i?mag values for sources belonging to an hypothetical 
sample characterised by Si.4GHz ^ 1 A*Jy (Sect. [5~2.2|l . 

The three methods give consistent results, which is 
reassuring and shows that we have a reasonable handle on 
the magnitudes of faint radio sources. The scaled R m a g 
for radio-quiet AGN are fainter than those derived for 
"typical" sources and are beyond the -R = 1.4 line. This 
is because these sources have the highest redshifts in the 
VLA-CDFS sample, and therefore their observed magni- 
tudes are more affected by K-correction effects. 

The horizontal dot-dashed line at -R ma g ~ 29.3 rep- 
resents the limit of the dee pest optical data cur rently 
available, the Hubble UDF (jBeckwith et al.ll2006l ). cov- 
ering about 11 arcmin 2 . The other horizontal lines indi- 
cate, from top to bottom, the limiting magnitudes for the 
Panoramic Survey Telescope and Rapid Response System 
(PAN-STARRSJ3, which from Hawaii will survey about 
3/4 of the sky down to R m a g ~ 26 during 10 years of oper- 
ation and the Large Synoptic Survey Telescope (LSST0 
which will be located in Chile and will provide a survey of 
about half the sky down to R ma g ~ 27.5 during 10 years 



13 http:/ /pan-starrs. ifa.hawaii.edu/ 

14 http://www.lsst.org 
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Figure 3. ifmag vs. the 1.4 GHz radio flux density for faint radio sources. Diagonal lines represent different values of R = 
log(S'i.4GHz/5'R, maB ), ranging from —1 (top) to 4 (bottom). The diagonal dashed line at R = 1.4 indicates the maximum value 
for SFG and the approximate dividing line between radio-loud and radio-quiet AGN, with SFG and radio-quiet AGN expected 
to populate the top left part of the diagram. The scaled R magnitudes of prototypical representatives of the three classes at 
•S'l.4GHz = 1 A*Jy are a l so shown, with SFG split into starbursts, spirals, and dwarfs. Finally, the mean radio and -f? mag values for 
sources belonging to an hypothetical sample characterised by 5i.4GHz ^ 1 MJy, as extrapolated from the VLA-CDFS sample, with 
error bars indicating the standard deviation, are also marked. The horizontal dot-dashed lines indicate the approximate point- 
source limits of (from top to bottom): PAN-STARRS, LSST, JWST, the deepest Hubble Space Telescope surveys, and the E-ELT. 
Survey areas (for LSST and PAN-STARSS), S/N ratios and exposure times (for the E-ELT and JWST) and operation/launch 
dates for future missions, or best guesses at the time of writing, are also shown. See text for more details. 



of operation. Further below there are the limiting magni- 
tudes for point sources for the James Webb Space Tele- 
scope (JWSTj3 and for the European Extremely Large 
TelescopS (E-ELT0 for the given signal-to- noise (S/N) 



15 http://www.stsci.edu/jwst/ 

16 http: //www. eso.org/sci/facilities/eelt/ 

17 I am using the E-ELT as an example because it has the 



ratios and exposure times I. Both these telescope will op- 
erate in observatory mode. It is important to keep in mind 

largest mirror (42 m in diameter) amongst the three very large 
telescopes being planned. The other two are the Thirty Meter 
Telescope (TMT; [http://www.tmt.org/ll and th e Giant Mag- 
ellan Telescope (GMT; http://www.gmto.org/). 
18 The JWST limit comes from the JWST Web pages. The 
E-ELT limit was derived from the E-ELT Exposure Time Cal- 
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that all of these limits are only approximate, as none of 
these telescopes are in operation yet. 

Given the progression towards lower R values going 
from (non-dwarf) starbursts, to spirals, and to BCDs and 
dwarf irregulars (Sect. I3~l2. 1(1 . the typical magnitudes of 
faint radio sources at a given flux density will depend 
critically on which of these SFG sub-classes dominates 
the counts. 

The contribution of "normal" spirals to faint radi o 
number counts is not well known. iHopkins et al.l (|2000f l. 
by converting the optical LF to a radio LF adopting a 
constant ratio Piagh?./Lb = 1/3, predict that normal 
galaxies will outnumber starbursts at Si.4GHz ~ 1 A'Jy. 
Since J? is a function of radio power, this assumption is 
too simplistic. I have used a different approach based on 
the SFR division introduc ed above (Sec t. r3.2.1|) . Namely, 
I have split the SFG LF of lSadler etafl (|2002h into two at 
A.4GHz = 8.4 x 10 21 W Hz -1 and evaluated the n umber 
counts assuming the evolution of IHopkins! {2004). The 
result is that "normal" spirals should be more numerous 
than starbursts for Si.4GHz < 5 /iJy. 

Are dwarf galaxies also going to be relevant at the 
flux densities of interest here (> 10 nano Jy)? Although 
we do not have any information on their radio LF, one 
can ma ke an educa t ed gu ess by converting their opti- 
cal LF. iLeon et al.l (|2008T l find that radio and optical 
power for Sb - Sc galaxies extending down to Mb ~ — 16 
scale linearly; the scaling factor is consiste nt with the 
mean values for the dwarf galaxies studied bv lLerov et al.l 
(2005), which rea ch Mb ~ — 14. I have then transformed 
the dwarf LF of iBlanton et all (120051 ) to a radio LF, 
which covers the range 10 18 < Pi.4GHz < 3 x 10 20 W 
Hz -1 . Assuming the same evolution as SFG and in- 
tegrating as done above to z max = 3, one then finds 
that dwarfs outnumber non-dwarf SFG at Si.4ghz < 20 
nanojy. This result is quite robust towards extrapola- 
tions of the optical LF to fainter magnitudes (M x ~ —10) 
with the same slope (a — —1.52), while a steepening 
of the LF (a — —2.0) would lead to an increase in the 
flux density below which dwarfs dominate by a factor 
~ 3. Note also that any evolution smaller than the as- 
sumed one would push further down the contribution 
of dwarf galaxies at a given flux density. However, the 
range in radio power covered by this radio LF does not 
extend to the values re ached by some BCDs (Pi,4gh z > 
5 x 10 20 W Hz -1 : e.g.. iHunt. Bianchi fc Maiolind I2OO5I : 



Sar gsvan fc Wedmanlboogt ). I have then multiplied the 
scaling factor by 10 0,5 to take into account the some- 
what larger radio-to-optical flux density ratio of BCDs 
(see Sect. I3.2TT) . with the result that dwarfs should out- 
number non-dwarf SFG at 5i.4GHz < 0-3 ^Jy- Since not 
all dwarfs are BCDs, dwarf galaxies will become the most 
numerous constituents of the radio sky at flux densities 
in the likely range 20 nanoJy < Si.4GHz < 300 nanojy. 

Fig. [3] shows that, down to 1 /j,3y, most starburst- 
like SFG and radio-quiet AGN should be detected by 
the LSST, that is they will have a counterpart in a 
large area survey. To be more specific, based on the 



culator (ETC) using the Laser- Tomography/Multi-Conjugate 
Adaptive Optics option. 



VLA-CDFS extrapolation, ~ 66% of starburst galax- 
ies with Si.4GHz zZ 1 /^Jy should have P mag < 27.5. 
This number represents an upper limit because by scal- 
ing the VLA-CDFS magnitudes I have assumed that the 
mean redshift is unchanged, while IPadovani et al.l (|2009t ) 
have found a strong correlation between redshift and 
magnitude. Therefore, K-correction effects will be larger 
and observed magnitudes will be fainter. On the other 
hand, "normal" spirals should outnumber starbursts for 
Si.4GHz < 5 fiJy. If even distant spirals are characterised 
by rest-frame R < 0.1 — 0.3, as is the case for local 
sources (Sect. I3.2.1|l . then pijy sources could be « 2 
magnitudes brighter than expected in the case of star- 
bursts and therefore within reach of PAN-STARRS and 
certainly of the LSST. For fainter radio samples optical 
magnitudes should get fainter, unless dwarf galaxies, with 
their low R values, become dominant, which would lead 
yet to another "brightening". For example, for R ~ —0.5 
and 5i.4GHz ~ 100 nanojy, P mag « 25. 

FR Is have on average larger R values but also higher 
radio flux densities, so their mean magnitude in an hy- 
pothetical 5i.4GHz 1 A*Jy sample is brighter than those 
for "typical" sources with 5*i.4ghz = 1 mJ v j although the 
R values are similar. Low-power ellipticals, however, will 
have smaller R and therefore brighter magnitudes. 

Sources having P mag > 27.5 will obviously be within 
reach of JWST and ELTs, which however will be covering 
a relatively small field of view (up to a few arcmin 2 ). 
Depending on the actual relative fraction of starbursts, 
spirals, and dwarfs, these observatories could be the main 
(only?) facilities to secure optical counterparts of nanojy 
radio sources. 

To address the question of what fraction of sources in 
the LSST and PAN-STARRS surveys will have an SKA 
counterpart requires a knowledge of the radio proper- 
ties of very faint optical sources, which at present we do 
not possess. Nevertheless, one can make some educated 
guesse s. The UDF number counts of Beckwi th et al.l 
(|2006T l show that for zab < 26 the surface density is 
~ 3 x 10 5 deg -2 , which implies, based on Fig. [T] that 
the majority of the objects beyond this limit are star- 
forming systems. SFG in the VLA-CDFS sample have 
observed R w 1.3, which means Si.4GHz ~ 3 ^tJy and 
« 70 nanoJy at P mag ~ 26 and ~ 27.5 respectively, that 
is within reach of the SKA. Optical selection however 
will be biased towards smaller R and therefore fainter ra- 
dio flux densities. More importantly, spirals and dwarfs, 
which might be the most numerous sub-classes, with their 
lower R values will have fainter radio flux densities for a 
given magnitude. As regards the small minority of radio- 
quiet AG N, they can reach R ~ —1.5 in the AGN cat- 
alogue of IPadovani et all l|l997l ). which would imply ex- 
tremely l ow radio flux densit ies even for PAN-STARSS. 
However. Iwhite et all l|2007h . by stacking FIRST radio 
images for ~ 40, 000 SDSS quasars, have shown that R 
increases with magnitude, which would imply higher ra- 
dio flux densities at faint magnitudes. In summary, the 
bulk of LSST and PAN-STARRS sources might not have 
radio flux densities within reach of the SKA, but at this 
point in time one cannot be more specific. 
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3.3 Mid-Infrared Band 

3.3.1 Typical mid-infrared-to-radio flux density ratios 

The infrared and radio emission are strongly and linearly 
correlated in SFG, d efining what is kno wn as the "IR- 
radio relation" (e. g. JSargent et al.1l2010l . and references 
therein). The most recent determination of the so-called 
(K-corrected) ?24 parameter, based on the COSMOS field 
and using an IR/radio-selected sample to reduce selec- 
tion biases, is q>24 = log^^m/Si^GHz) = 1-26 ± 0.13 
l|Sargent et alJuOld ). As mentioned above (Sect. [3T.1[I . 
radio selection will favour objects with larger radio flux 
densities, and therefore smaller 524, but given the tight- 
ness of the r elation the de c rease is only ~ 0.3 (based 
on Tab. 3 of ISargent et alj l|2010l )'), BCDs also appear 
to ha ve q parameters globally consis tent with tho se of 
SFG dHunt. Bianchi fe Maiolindl2005h . even though iBelll 
l|2003h argues that both the IR and radio luminosities 
of dwarf galaxies significantly underestimate the SFR by 
similar amounts. Radio-quiet AGN have 524 values sim- 
ilar to those of SFG, while radio-loud ones have much 
smaller ones. For this band we cannot extrapolate from 
the VLA-CDFS sample since most of the sources are un- 
detected in the 24 /im band in publicly available cata- 
logues. 

Evolution in the mid-infrared and radio band is 
broadly similar for SFG (Ranalli, Comastri & Sctti 2005) 
and AGN (e.g., IWeedman fc Houckl I2009D . Therefore, 
typical flux density ratios will be largely unaffected by it. 
No information is available on the mid-infrared evolution 
of FR Is, which in any case are very weak IR emitters. 



3.3.2 Mid- Infrared flux densities of faint radio sources 

Figure U plots the 24pim vs. 1.4 GHz radio flux density 
and shows the expected /24 M m for "typical" radio-quiet 
AGN and SFG (from NED), scaled to 1 ^Jy (the values 
for two "typical" dwarf galaxies are fully consistent with 
these and are not included for clarity) and the locus of 
SFG and radio-quiet AGN (Sect. [3X1]) • The two methods 
give perfectly consistent results. The two "typical" FR Is 
have /24 M m ~ 0.2 — 0.8 nanojy and are therefore way off 
the plot. 

The two top horizontal dot-dashed lines represent 
the point-source limits for the Wide-field Infrared Survey 
Explorer (WISEj3 and AKAR0 

Both are performing 
all-sky infrared surveys. However, I give the AKARI limit 
for pointing mode as the survey limit for the band clos- 
est to the one of interest (18 (im) is off the plot (120 
mjy). The horizontal line at /24 Mm ~ 40 /iJy indicates 
the deepest mid-infrared data currently available, that 
is the Spitzer observat ions of the GOODS/FIDEL field 
(|Bethermin et al.ll2010T ), covering about 0.23 deg 2 . The 
bottom two horizontal lines denote the limiting flux den- 
sities for point sources for JWST and the SPace Infrared 
telescope for Cosmology and Astrophysics (SPICA j3 f° r 



http:/ /wise. ssl.berkeley.edu/ 

http: / /www. ir.isas.jaxa.jp/ ASTRO-F / 

http: / /www. ir.isas.jaxa.jp/SPICA/ 



the given S/N ratios and exposure times. Both these tele- 
scope will operate in observatory mode. These two limits 
are only approximate, as none of these telescopes are in 
operation yet. 

Figure |3] shows that most SFG and radio-quiet AGN 
should be detected by JWST down to ~ 1 pjy and by 
SPICA down to ~ 100 nanojy. Both telescopes will be 
covering a relatively small field of view (up to ~ 15 
arcmin 2 ). Only a tiny fraction of /iJy and nanojy SFG 
and radio quiet AGN will have a counterpart in an all-sky 
infrared surveys. As regards FR Is, they will be largely 
undetected even by SPICA. 

On the other hand, all extragalactic sources in 
the WISE and AKARI all-sky surveys will easily have 
an SKA counterpart, as they will be characterised by 
SlAGHz > 100 /LtJy. 



4 REDSHIFTS OF FAINT RADIO SOURCES 

A vital component in the identification of astronomical 
sources is redshift, which allows powers to be estimated 
and LFs to be derived. Figure [S] plots i? mag vs. 1.4 GHz 
radio flux density but, at variance from Fig. [3] the marked 
limits refer to spectroscopic or photometric redshifts. 

The top horizontal dot-dashed line represents the 
current limiting magnitude for "classical" photometric 
redshifts, which rea ch -R mag ~ 24 and z » 1.5 (e.g., 
iMainieri et al.1 [2008); photometric redshifts for special 
classes of sources, e.g., z ~ 7 drop -outs can be derive d 
with HST down to Fab ~ 28 (jBouwens et al.l 120101 ). 
Long exposures (~ lOh) with 8/10 m telescopes can se- 
cure spectroscopic redshifts in the case of strong emission 
lines down to 7? m ag ~ 26, so this value represents a hard 
limit. Finally, indicative limiting magnitudes for obtain- 
ing spectra for point sources with the E-ELT and JWST 
are also shown 22 !. 

Fig. shows that, down to 1 fj,Jy, only approxi- 
mately half of the starburst-like SFG will be within reach 
of 8/10 m telescopes in terms of obtaining a spectro- 
scopic redshift. To be more specific, based on the VLA- 
CDFS extrapolation, ~ 50% of starburst galaxies with 
5i.4GHz > 1 A'Jy should have -R mag > 26 and ~ 60% will 
be above -R mag = 25, probably a more realistic limit for 
a spectrum. As discussed above ("Sect. [3"T2.3[1 . these frac- 
tions are robust upper limits. On the other hand, spiral, 
and maybe dwarf galaxies, should be more common than 
starbursts at fijy and nano-Jy levels, respectively, which 
would lead to brighter magnitudes, thereby making the 
job of deriving a redshift easier. Future facilities like the 
E-ELT and JWST will obviously allow the determina- 
tion of spectroscopic redshifts for fainter (i? mag > 26) 

22 The JWST limit refers to a resolution of 1,000 and A = 
1 £tm and comes from the JWST Web pages. The E-ELT lim- 
its apply to a resolution of 1,000, the R band, and was derived 
from the E-ELT ETC using the Laser- Tomography/Multi- 
Conjugate Adaptive Optics option. Since the choice of instru- 
ments for the E-ELT has not been finalised yet, these limits are 
very preliminary. Note that in both cases the magnitudes refer 
to the continuum: the presence of emission lines will obviously 
lead to shorter exposure times. 
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Figure 4. 24/im flux density vs. the 1.4 GHz radio flux density for faint radio sources. The diagonal dashed line represents the 
locus of SFG and radio-quiet AGN based on the "IR-radio relation" . The scaled IR flux densities of prototypical representatives 
of the three classes at Si.4GHz = 1 ^Jy are a l so shown, with FR Is being so faint as to be actually off the plot at /24 Mm ~ 0.2 — 0.8 
nanojy. The horizontal dot-dashed lines indicate the approximate point-source limits of (from top to bottom): WISE, AKARI 
(pointing mode), the deepest Spitzer surveys, JWST, and SPICA. Launch dates for future missions, or best guesses at the time of 
writing, are also shown. See text for more details. 



objects, albeit with relatively long exposures. It might 
be also reasonable to suppose that the LSST will reach 
the required accuracy to determine photometric redshifts 
for sources one magnitude brighter than its 10 yr sur- 
vey limit, thereby reaching i? mag ~ 26.5. However, both 
LSST and PAN-STARRS will work with photometric 
systems with ~ 6 optical broadband filters, similar to 
those traditionally used in astronomy, which means that 
they will be more prone to color-redshift degeneracies 
than systems using a larger number of narrower filters 
(|Benftez et alj|2009l ). Since the reddest filter is going to 



be in the y-band (~ 1.05 /im), this will also limit the 
maximum achievable redshift. 

Note that the magnitude limits of the current deep- 
est large-area spectroscopic surveys are r ^ 17.8 (galax- 
ies) and i < 19.1 (quasars) for the SDSS^and bj < 20.85 
for the 2dF QSO redshift survejH 

In summary, many of the 5*i.4ghz 1? 1 sources 
might be too faint for 8/10 m telescopes to be able to 
provide a redshift and the situation might get worse at 

23 http://www.sdss.org/dr7/ 

24 http://www.2dfquasar.org/ 
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Figure 5. ifmag vs. the 1.4 GHz radio flux density for faint radio sources. Diagonal lines represent different values of R = 
log(S'i.4GHz/5'R, maB ), ranging from —1 (top) to 4 (bottom). The diagonal dashed line at R = 1.4 indicates the maximum value 
for SFG and the approximate dividing line between radio-loud and radio-quiet AGN, with SFG and radio-quiet AGN expected 
to populate the top left part of the diagram. The scaled R magnitudes of prototypical representatives of the three classes at 
•S'l.4GHz = 1 A*Jy are a l so shown, with SFG split into starbursts, spirals, and dwarfs. Finally, the mean radio and i? mag values for 
sources belonging to an hypothetical sample characterised by Si.4GHz ^ 1 A*Jy, as extrapolated from the VLA-CDFS sample, are 
also marked. The horizontal dot-dashed lines indicate the approximate limiting magnitudes to derive photometric and spectroscopic 
redshifts. From top to bottom: the current limit for photometric redshifts, the JWST and an E-ELT limit (for the given S/N ratios 
and exposure times), the faintest magnitude for which a spectroscopic rcdshift can be obtained with 8/10 m telescopes, a possible 
limit for photometric redshifts with the LSST, a fainter E-ELT limit, and the limit of photometric redshifts for z ~ 7 drop-outs 
with the Wide Field Camera 3 (WFC3) on board HST. See text for more details. 



fainter flux densities, unless dwarf galaxies take over. This 
means that JWST and the ELTs might be the main fa- 
cilities to secure redshifts of /^Jy radio sources. But even 
they could have problems in the nanojy regime. 



5 DISCUSSION 

The simple approach employed in this paper suggests 
that (non-dwarf) SFG play a very important role in the 
sub-/ijy sky, something which had been already realised 
by many authors. What is new here is that also low- 
power (Pi.4GHz < 10 20 W Hz" 1 ) ellipticals but espe- 
cially dwarf galaxies are very relevant, with the latter 
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very likely being the most important constituent of the 
faint radio sky. The arguments I have used to reach this 
conclusions appear quite robust. Indeed, by limiting the 
integration to z max = 3 for lack of information, the de- 
rived values on the surface densitie s of these classes are 
lower limits. IChristlein et al.l ((2009) have recently used a 
new maximum likelihood method to derive a LF down to 
Mr ~ —14, which disp lays a faint end slop e a = —1.62 
steeper than the LF of iBlanton et alj (|2005l ) I have used 
for my calculations (which has a = —1.52). The number 
(and therefore surface) density of dwarf galaxies is then 
probably even higher than I have assumed. 

The radio, far-infrared, and X-ray bands all trace 
star formation, which implies relatively tight relation- 
ships between these three bands. This means that we can 
reasonably predict the mid-infrared and X-ray proper- 
ties of SFG with faint radio flux densities. The situation 
is different in the optical band, which is dominated by 
emission from evolved stars. Since "normal" spirals and 
dwarf galaxies have radio-to-optical flux density ratios 
smaller than starbursts because of their lower star for- 
mation rate, they will be brighter in the optical band at 
a given radio flux density ratio. The optical properties 
of faint radio sources will then depend on the population 
mix of starbursts, spirals, and dwarf galaxies, which at 
present cannot be determined very accurately. Given our 
ignorance of the radio LF and evolution of dwarf galax- 
ies, in fact, it is not straightforward to predict at which 
flux densities they will become the dominant population. 
By converting the optical LF to a radio LF, under the 
assumption (not inconsistent with the available data) of 
a linear proportionality between the two powers and of 
an evolution similar to that of SFG, dwarfs should out- 
number non-dwarf SFG in the 20 - 300 nanojy range. 
However, this cannot be the whole story, as the esti- 
mated radio LF reaches only Pi.4GHz ~ 10 21 W Hz -1 , 
while there are BCDs with P i. 4 ghz up to w 10 22 W Hz -1 
(jSargsvan fc Wedmarj2009h . There has then to be a high- 
power component of the LF, which will translate into a 
higher number of sources at high flux densities, which at 
present cannot be quantified. 

Can the two "new" populations be re levant also 
for the radio background? ISingal et al.l (|2010l ) have con- 
cluded that the main background contributors have to be 
faint (sub-pjy) radio sources and suggested these to be 
ordinary star-forming galaxies at z > 1 characterised by 
an evolving radio far-infrared correlation, which increases 
towar d the radio loud with redshift (but see lSargent et all 
(2010) for evidence that the local radio far-infrared rela- 
tion still holds up to z ~ 5). iMassardi et alj l|2010l ) (see 
also Padovani et al., in preparation) have shown that pop- 
ulation synthesis models of the radio sky, which currently 
involve only "classical" radio sources, fail to explain by 
a factor of a few (see th eir Fig. 9) the very recent back- 
ground measurements of lFixsen et al . (2010). Since back- 
ground emission from a population of objects is equal 
to J SN(S)dS (where N(S) are the differential counts), 
which is a combination of surface and flux density, the 
right mix of numbers and fluxes are needed to provide 
a relevant contribution. Therefore, while low-power ellip- 
ticals, which are both faint and not that numerous, are 
most likely unimportant, for dwarf galaxies the situation 



is less obvious. On one hand, even with the largest scal- 
ing factor used in Sect. 13.2.31 dwarf galaxies account for 
only « 8% of the sub-mjy background. But on the other 
hand, the radio LF is not known and, as discussed above, 
the radio number counts of dwarf galaxies could well be 
higher. 

The determination of the surface densities for all 
classes of radio sources allows us also to determine the 
intrinsic ratio between radio-loud and radio-quiet AGN. 
The surface density of SSRQs and FSRQs combined is 
~ 6.3 deg -2 , while that of radio-quiet, broad-lined AGN 
is Ri 5 x 10 3 deg" 2 , since unabsor bed AGN are sup posed 
to make up ~ 1/5 of the total (jGilli et alJl2007t >. This 
translates to a radio-loud fraction in broad-lined AGN 
~ 0.1%. If one considers both broad- and narrow-lined 
sources, which means including also FR lis in the radio- 
loud class and all radio-quiet AGN, one gets a similar 
fraction « 0.2%. These numbers are much smaller than 
the oft-quoted value of ~ 10 — 20%. However, it has been 
known for some time that the rad io-loud fraction drops 
with dec reasing optical lum inosity (|Padovanilll993l ). Very 
recently, Ijiang et alj (|2007h have shown that the fraction 
of radio-loud quasars at z — 0.5 declines from 24% to 6% 
as luminosity decreases from M2500 = —26 to —22, while 
at M2500 = — 26 this fractions declines from 24% to 4% as 
redshift increases from 0.5 to 3. Both these results suggest 
that the intrinsic, global ratio has to be < 5%. But if the 
radio-loud fraction were this large, the radio-loud quasar 
surface density should be ~ 250 deg -2 , that is ~ 40 times 
larger than estimated here and also ~ 5 larger than that 
of FR lis, which are supposed to be the parent popula- 
tion of both SSRQs and FSRQs. Scaling up the FR II 
surface density as well to avoid this paradox, one would 
end up having more FR lis than FR Is, which would not 
make any sense at all, since FR lis are more powerful. In 
short, the surface densities derived in this paper point to 
an AGN radio-loud fraction much smaller than normally 
derived. This difference is likely to be ascribed to two 
factors: a) the usually quoted value has been obtained 
for b right, optically selected samples (|Kellermann et all 
1989), which include radio-loud quasa rs with their optical 
flux boosted by relativistic beaming |Goldschmidt et alj 
1999), which artificially increases their fraction; b) radio- 
loud AGN are on average more powerful than radio-quiet 
ones (e.g.. IZamfir. Sulentic fc Marz iani 2008). When one 
reaches the very faint end of the optical LF only radio- 
quiet AGN will be present and therefore the integrated 
radio-loud fraction will be quite small. In other words, the 
usually quoted value refers to the bright part of the LF 
and, when integrated over the full range of powers, the 
resulting radio-loud fraction is much smaller. The second 
factor is likely to be the most relevant one. 

The issue of source confusion goes beyond the scope 
of thi s paper. I will just point out that IWindhorst et all 
(|2008l ) have shown that ultra-deep radio and optical sur- 
veys may slowly approach the natural confusion limit, 
where objects start to overlap with their neighbours due 
to their finite sizes and not because of the finite instru- 
mental resolution, which causes the instrumental confu- 
sion limit. Based on the model of lWindhorstl (I2003T ). con- 
fusion will not be a problem for the SKA if source sizes are 
< 1.6" at Siagtu = 1 ^Jy (or < 0.25" at Si. 4 ghz = 10 
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nano Jy) and the instrumen tal resolution is commensu- 
rate (Windhorst et al.l[20ol ). However, since IWindhorstj 
l|2003 ) did not include radio-quiet AGN, low-power ellip- 
ticals, or dwarf galaxies, confusion l imits are bound to be 
higher. Indeed. lOwen fc Morrison] l|2008l ). based on the 
relatively flat number counts derived from their deep ra- 
dio data, have suggested that the natural confusion limit 
may be reached near 1 /ijy. 



6 CONCLUSIONS 

I have used a simple method to work out the source pop- 
ulation of the sub-/ijy sky, which relies on estimating the 
smallest flux density and the largest surface density for 
all classes of radio sources. These two parameters can 
be derived quite robustly from the local minimum ra- 
dio power and number of sources per unit volume, max- 
imum redshift, and luminosity and density evolution. I 
have then estimated for the first time the X-ray, optical, 
and mid-infrared fluxes these faint radio sources are likely 
to have by using prototypical sources, typical flux density 
ratios, and extrapolations from the VLA-CDFS sample. 
Prognosticating these multi-wavelength properties is ex- 
tremely relevant for the identification of /xJy and nanoJy 
radio sources and to maximise the synergy between the 
SKA and its pathfinders with future missions in bands 
other than the radio. My main results can be summarised 
as follows: 

(1) The sub-/ijy sky should consist of radio-quiet AGN 
(~ 2 x fO 4 deg -2 ) and star-forming galaxies (« 3 x 10 6 
deg -2 ), both of which should get to Si.4GHz ~ 0.1 — 1 
nanoJy. In agreement with previous studies, I find that 
classical, powerful radio sources, that is radio quasars and 
FR lis, do not make it to sub-/^Jy flux densities and reach 
5i.4GHz ~ 0.1 mjy, with BL Lacs and FR Is getting to 
the pa 1 — 10 /iJy level. 

(2) The intrinsic fraction of radio-loud AGN, inte- 
grated over the whole bolometric luminosity function, 
is « 0.1 — 0.2%, that is about two orders of magnitude 
smaller than the oft-quoted value of « 10 — 20%. The 
latter value refers only to the bright part of the opti- 
cal luminosity function and is also biased because of the 
likely flux boosting in radio-loud quasars due to relativis- 
tic beaming. 

(3) Two "new" populations, which have not been con- 
sidered previously, appear to be very relevant: low-power 
(fi.4GHz < 10 20 W Hz -1 ) ellipticals and dwarf galax- 
ies. Using the available (scanty) information, the former 
should reach similar flux and surface densities as radio- 
quiet AGN, while the latter could easily be the most 
numerous component of the faint radio sky (> 5 x 10 6 
deg -2 ), with flux densities as low as ~ 1 nanoJy. Since 
most galaxies at the faint end of the LF are blue, most 
dwarfs in the Universe (and therefore most radio sources) 
should be of the star-forming type. While low-power el- 
lipticals are most likely unimportant contributors to the 
radio background, the verdict is still open for dwarf galax- 
ies. This issue is quite important also in light of the very 
recent background me asurements reporte d by the AR- 
CADE 2 collaboration (jFixsen et alj|2010l ). 



(4) The bulk of the /^Jy population, which is then most 
probably going to be made up of star-forming galaxies, 
is likely to have X-ray fluxes beyond the reach of all cur- 
rently planned X-ray missions, including IXO. The same 
applies to FR Is. Even IXO will only detect radio-quiet 
AGN with radio flux densities > 1 /zJy. The situation 
will obviously be worse in the nanoJy regime, where very 
few radio sources will have an X-ray counterpart in the 
foreseeable future. On the other hand, basically all extra- 
galactic sources in the eRosita All-Sky and Wide Surveys 
and WFXT Wide Survey will have an SKA counterpart, 
as they will be characterised by 5i.4GHz > 1 H-3y. This 
will help in the identification of, for example, the 10 mil- 
lion or so point sources expected in the latter, by also 
providing very accurate positions. 

(5) In the mid-infrared most star-forming galaxies and 
radio-quiet AGN should be detected by JWST down to 
radio flux densities ~ 1 pjy and by SPICA down to ~ 
100 nanoJy, while FR Is will be largely undetected even 
by SPICA. Both telescopes will however be covering a 
relatively small field of view. Only a minute fraction of 
fijy and nanoJy star-forming galaxies and radio quiet 
AGN will have a counterpart in current all-sky infrared 
surveys (WISE and AKARI). On the positive side, all 
extragalactic sources in the WISE and AKARI all-sky 
surveys will easily have an SKA counterpart, as they will 
be characterised by £i.4GHz > 100 /iJy. 

(6) Since the radio, far-infrared, and X-ray emission 
all trace star formation, this implies relatively tight re- 
lationships between these three bands, which in turn 
means that we can reasonably predict the mid-infrared 
and X-ray properties of star-forming galaxies with faint 
radio flux densities. The situation is different in the op- 
tical band, where evolved stars dominate. Objects char- 
acterised by higher star formation rates will have larger 
radio-to-optical flux density ratios, which means fainter 
magnitudes for a given radio flux density. The typical 
magnitudes of the optical counterparts of faint radio 
sources depend then on which type of star-forming galaxy 
(starburst, spiral, or dwarf) will be predominant. More- 
over, in the optical band K-correction and intergalactic 
absorption effects are important and will result in fainter 
than expected magnitudes. Assuming a maximum value 
of star formation rate for "normal" spirals, which trans- 
lates to a maximum radio power, I estimate that these 
sources should outnumber starbursts for Siaghz 5 /xJy. 
If even distant spirals are characterised by relatively low 
rest-frame radio-to-optical flux density ratios, then most 
/iJy sources should be detected by PAN-STARRS and 
certainly by the LSST. At fainter radio flux densities 
optical magnitudes should also get fainter, unless dwarf 
galaxies, with their lower radio-to-optical flux density ra- 
tios, become dominant. Depending on the relative frac- 
tion of starbursts, spirals, and dwarfs, JWST and espe- 
cially the ELTs could be the main, or perhaps even the 
only, facilities capable of securing optical counterparts of 
nanoJy radio sources. On the other hand, and for the 
same reasons discussed above, the bulk of PAN-STARRS 
and LSST sources might not have radio flux densities 
within reach of the SKA. 

(7) As regards redshifts, the same complications de- 
scribed above apply, since the optical band is still in- 
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volved. Within the same uncertainties, then, many of the 
sources with Si.4GHz 1? 1 fiJy might be too faint for 8/10 
m telescopes to be able to provide a redshift determi- 
nation and the situation might get worse at fainter flux 
densities, unless dwarf galaxies are prevalent. This means 
that JWST and particularly the ELTs might be the pri- 
mary facilities to secure redshifts of fiJy radio sources. 
But even they could have problems in the nanoJy regime. 

In summary, the SKA and its pathfinders will have a 
huge impact on a number of open problems in extragalac- 
tic astronomy. Apart from the "obvious" study of star- 
forming galaxies in their various incarnations and "clas- 
sical" radio sources, these include also less evident ones, 
ranging from what makes a galaxy radio-loud (through 
the study of low-power ellipticals), to why most AGN 
are radio-quiet (by selecting large samples independent 
of obscuration), to the incidence and evolution of dwarf 
galaxies (by providing a "cleaner" radio selection, which 
might by-pass the surface brightness problems of optical 
samples), to the resolution of the radio background. 

Identifying faint radio sources, however, will not be 
easy. On large areas of the sky the SKA will be quite alone 
in the multi-wavelength arena, with the likely exception 
of the optical band and even there probably only down 
to « 1 ^Jy. SPICA, JWST, and especially the ELTs will 
be a match for the SKA but only on small areas and 
above 0.1 — 1 /iJy. At fainter flux densities one might 
have to resort to "radio only" information, that is HI 
redshifts, size, morphology, spectral index, etc., although 
I think this will not be sufficient. On the bright side, most 
sources from currently planned all-sky surveys, with the 
likely exception of the optical ones, will have an SKA 
counterpart. 
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